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Abstract 

We solve the Faddeev equation in an exactly Poincare invariant formulation of the three-nucleon problem. 
The dynamical input is a relativistic nucleon-nuclcon (NN) interaction that is exactly on-shell equivalent to 
the high precision CD Bonn NN interaction. S'-matrix cluster properties dictate how the two-body dynamics 
is embedded in the three-nucleon mass operator (rest Hamiltonian) . We find that for neutron laboratory 
energies above « 20 MeV relativistic effects on Ay are negligible. For energies below « 20 MeV dynamical 
effects lower the nucleon analyzing power maximum slightly by ~ 2% and Wigner rotations lower it further 
up to w f 0% increasing thus disagreement between data and theory. This indicates that three-nucleon forces 
(3NF) must provide an even larger increase of the Ay maximum than expected up to now. 

PACS numbers: 21A5.+V, 24.70. +s, 25.10.+S, 25.40.Lw 
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I. INTRODUCTION 

iigh precision nucleon-nucleon (NN) interactions such as AV18 CDBonn Q], Nijm I, II and 
93 ^ accurately describe the NN data set up to about 350 MeV. When these interactions are used 
to predict binding energies of three-nucleon (3N) systems they underestimate the experimental 
bindings of and '^He by about 0.5-1 MeV This missing binding energy can be cured by 

introducing a three-nucleon force (3NF) into the nuclear Hamiltonian Q]. 

The study of elastic nucleon-deuteron (Nd) scattering and nucleon-induced deuteron breakup 
also revealed a number of cases where the nonrelativistic description based on pairwise interactions 
is insufficient to explain the data. Generally, the studied discrepancies between a theory based 
on NN interactions only and experiment become larger with increasing energy of the 3N system. 
Adding a 3NF that includes long-range 27r exchange to the pairwise interactions leads in some cases 
to a better description of the data. The parameters of such a 3NF must be separately adjusted to 
the experimental binding of and ^He [3, Q, for each NN interaction. The elastic Nd angular 
distribution in the region of its minimum and at backward angles is the best studied example 6^\. 
The clear discrepancy in these angular regions at energies below ~ 100 MeV nucleon laboratory 
energy between a theory based on NN interactions only and the cross-section data can be removed 
by adding modern 3NFs based on chiral effective field theory [ol to the nuclear Hamiltonian. At 
energies higher than ~ 100 MeV current 3NFs l3, U] only partially improve the description 
of cross section data and the remaining discrepancies, which increase with energy, indicate the 



possibility of relativistic effects 
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14l | . The need for a relativistic description of 3N scattering 
was also raised when precise measurements of the total cross section for neutron-deuteron (nd) 
scattering [l^ were analyzed within the framework of nonrelativistic Faddeev calculations 16]. 
NN interactions alone were insufficient to describe the data above ~ 100 MeV. 

In few-body models off-shell effects, relativistic effects, and three-body force contributions can- 
not be cleanly separated. This is because different two-body interactions that give the same two- 
body S matrix „,a.e. by a unitary scattering equivalence H^To maintain this equivalence 
at the three-body level requires additional three-body interactions [18] in one of the Hamiltonians. 
Since relativistic two-body models are fit to the same data as the corresponding nonrelativistic 
models, there is a similar on-shell two-body scattering equivalence. While the relativistic and non- 
relativistic three-body predictions will be different one can in principle make up the difference (in a 
chosen frame) with a suitable three-body interaction. So while it is possible to simulate "relativistic 
effects" with a three-body interaction, a Poincare invariant treatment of the dynamics provides the 
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most direct way to model the consequences of imposing Poincare invariance and S-matrix cluster 
properties. 

In this paper we investigate one particular representation of the Poincare invariant three-body 
problem. We compare the predictions of relativistic and nonrelativistic three-body calculations 
where the input two-body interactions give the same two-body S matrix, have the same inter- 
nal two-body wave functions, and have a kinematic three-dimensional Euclidean symmetry. The 
relativistic and nonrelativistic models differ in how these interactions appear in the three body- 
problem. In addition, the internal and single particle variables are related by Galilean boosts in 
the nonrelativistic case and Lorentz boosts in the relativistic case. In the nonrelativistic case the 
two-body interactions for each pair are added to the center of mass kinetic energy operator. In 
the relativistic case, the non-linear relation between the two and three-body mass operators must 
be respected in order to obtain a scattering matrix that clusters into a product of the identity and 
the input two-body 5-matrix [19| . This non-linear dependence of the three-body invariant mass 



operator on the two-body interaction has dynamical consequences for the three-body system which 
complicates the structure of the Faddeev kernel. 

In {20! we used a Poincare invariant formulation of the 3N scattering problem. A technique for 
constructing the relativistic nucleon-nucleon interaction from a standard high-precision interaction 

n 

was given in [21]. We used the same technique to construct the transition operators that appear 
in the kernel of the relativistic Faddeev equation. Application to a 3N bound state supported the 
relativistic effects previously found in 22 1. 

Realistic NN interactions are fit by properly transforming experimental data to the center of 
momentum frame and fitting 5-matrix elements computed using the nonrelativistic Schrodinger 
equation to this data. While the same data could be precisely fit using 5-matrix elements computed 
from a relativistic Schrodinger equation, this has not been done 
to construct realistic interactions. 



23] with the same precision used 



In [2l|], instead, an analytical scale transformation of momenta was used to relate NN interac- 
tions in the nonrelativistic and relativistic Schrodinger equations in such a way, that the two-body 
scattering matrix elements are identified as Snr{Ecm) = Sr{Ecm) = Sexp{Ecm) as functions of the 
center of momentum energy [24] . In this work we use an alternative procedure [3] that generates 
a relativistic nucleon-nucleon interaction with the property that the relativistic and nonrelativistic 
two-body 5 matrices satisfy 5r(k^) = 5nr(k^) = 5exp(k^), where 5exp(k^) is the experimental 
two-body 5 matrix, and k^ is the cm momentum of one of the particles. 

When high-precision potentials are determined [2^] by properly Lorentz transforming scattering 
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data from the laboratory frame to the center of momentum frame the Lorentz invariant scalar 
product ptarget ' Pbeam = rnEj, = 2k^ + is used to relate the laboratory beam energy Ej, to 
the cm. momentum k^. The potential is determined by comparing the transformed experimental 
scattering observables to the scattering observables computed using the nonrelativistic Lippmann- 
Schwinger equation, identifying the appearing in the Lippmann-Schwinger equation with the 
k^ computed from the invariant ptarget 'Pbeam- With this procedure the resulting interactions are 
constructed so the S-matrix elements in the relativistic and nonrelativistic cases are identified as 
functions of k^ rather than cm energy. Though the difference in the two approaches leads to a 
small mismatch in the relativistic and nonrelativistic momentum, the interactions generated by the 
analytic scale transformation provide a useful first step to investigate the effects of the non-linear 
relation between the three-b ody mass operator and two-body interactions. 



In our initial studies 



3, 
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28l | the interaction generated by the analytic scale transformation 



was used to study the changes in elastic nd scattering and breakup observables when the nonrela- 
tivistic form of the kinetic energy is replaced by the relativistic one and a proper treatment of the 
dynamics is included. We found that the elastic scattering cross section is only slightly influenced 
by relativity. Only at backward angles and higher energies the elastic cross sections are increased 
by relativity. 

Due to the selectivity of the breakup reaction, however, regions of phase-space were found at 
higher energies of the incoming nucleon where relativity leads to a characteristic pattern by which 
relativity changes the nonrelativistic breakup cross section. Namely, in this region of phase-space 
fixing the angle of the first detected nucleon and changing the angle of the second nucleon provides 
variations of the nonrelativistic cross section by relativity, increasing or decreasing it by a factor 
up to ~ 2. For spin observables the implemented relativistic features lead only to small effects. 

Recently an interesting estimate of "relativistic corrections" has been performed and its effect 
on low-energy nd analyzing power Ay has been estimated using the plane-wave impulse approx- 



imation 29(1 . The estimate is based on a perturbative realization of the Poincare Lie algebra to 
leading order in l/cP [30]. The calculations are done in the plane- wave impulse approximation. A 
large increase by ~ 10% of the Ay maximum at laboratory energy En = 3 MeV has been found. 
The authors comment that their estimates are both exploratory and incomplete. In addition to the 
absence of final state interaction, there are a number of other important differences with an exact 
formulation of this problem. Such a large effect, which would significantly reduce the discrepancy 
between theory and data in the region of the Ay maximum, calls for a relativistic study in an 
exactly Poincare invariant treatment of the three nucleon dynamics. 
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Our previous study [20[], performed without inclusion of Wigner rotations, is too limited for 
spin observables. Therefore, in order to make definite conclusions for Ay we perform a complete 
dynamical calculation including the effects of Wigner rotations. We focus on that issue and do not 
include 3NF's. 

The paper is organized as follows. In section II we discuss the construction of our Poincare 
invariant dynamical model. This includes a discussion of how high-precision interactions are used 
to construct the three-body mass operator (rest Hamiltonian). In section III we discuss spin 
observables in Poincare invariant quantum theory. In Sec. IV we discuss the formulation of 
the Faddeev equation to construct scattering observables for this three-nucleon mass operator. 
This includes an exact treatment of the Faddeev kernel which avoids the approximations used 



m 



20|, 



27 



28( 1 . We solve the relativistic 3N Faddeev equation with and without Wigner spin 
rotations. We show and discuss results for the neutron analyzing power Ay. Sec. V contains a 
summary and conclusions. 



II. POINCARE INVARIANT DYNAMICS 



In a quantum theory the principle of special relativity requires that the probabilities computed 
for equivalent experiments done in different inertial coordinate systems are identical. Since inertial 
coordinate systems are related by Poincare transformations, it follows [31] that equivalent states in 
different inertial coordinate systems are related by a unitary representation, ^7(A, a), of the Poincare 
group. This emphasis on the invariance of experimental measurements in different inertial frames 
is different than the covariance requirements that are historically motivated by the way symmetries 
are realized in classical wave equations. 

Since any representation of the Poincare group can be decomposed into a direct integral of 
irreducible representations, one way to construct a Poincare invariant dynamics is to build it 
out of irreducible representations. The transformation properties of irreducible representations 
are well known and completely determined by group theoretical considerations. The dynamics is 
contained in the spectrum of the physical mass and spin operators which determines the values 
and multiplicities of the Casimir invariants that appear in this decomposition. 

Our construction begins with one-particle representations, which are irreducible representations. 
The particle's mass and spin fix the eigenvalues of the two Casimir invariants of the Poincare 
group. For computations it is necessary to choose a basis for the irreducible representation space. 
This is done by choosing a maximal set of commuting Hermitian functions of the infinitesimal 



5 



Poincare generators. In addition to the mass and spin, it is possible to find four additional mutually 
commuting non-invariant functions of the generators. There is a second set of four operators that 
are conjugate to the non-invariant commuting observables. These operators change the eigenvalues 
and determine the spectrum of the commuting observables. All ten generators can be expressed 
as functions of these eight non-invariant operators and the two Casimir operators. The irreducible 
representation space, is the space of square integrable functions of the eigenvalues of the four 
commuting operators [30]. 

Our choice of basis for irreducible representation spaces is the simultaneous eigenstates of the 
linear momentum p and the 3-component of the canonical spin, jcz, which is the observable corre- 
sponding to the spin measured in the particle's rest frame if the particle is transformed to its rest 
frame with a rotationless Lorentz transformation. In this basis the irreducible unitary representa- 
tion of the Poincare group is 



[/(A,a)|(j,m)p,/i) = m)p',/x')e'P'%/^I)^„,Ji?-Hp'MAi?(p/m)] (1) 



L0{p) 



where p' = Ap, uj{p) = \/p • p + rn?, and B^p/m) is the rotationless Lorentz transformation that 
takes a particle of mass m at rest to momentum p. The quantity R^{A,p) := B^^{p' /m)AB{p/m) 
is the standard rotationless-boost Wigner rotation. The representation ([1]) is unitary for states with 
a 6{p — p') normalization in the momentum. The important observation is that all mass m > 
spin j irreducible representations of the Poincare group in the {p,jcz} basis have this form. 

The two or three-nucleon Hilbert space is the tensor product of two or three single-nucleon 
irreducible representation spaces; 7i ®7i oi TL (^TL (^Ti. On each of these spaces 

UQ{A,a) = U{A,a)(^U{A,a) Uo{A,a) = U{A,a) U{A,a) U{A,a) (2) 

define kinematic representations of the Poincare group. These representations are reducible and 
do not contain any dynamics. We build dynamical irreducible representations by adding suitable 
interactions to the mass Casimir operator of non-interacting irreducible representations. The first 
step needed to introduce interactions is then to decompose these non-interacting tensor product 
representations into a direct integral of irreducible representations. This is accomplished with 
Poincare group Clebsch-Gordan coefficients in our chosen {p,jcz} basis. The Poincare group 
Clebsch-Gordan coefficients are the expansion coefficients of linear combination of tensor product 
states that transform irreducibly. The desired non-interacting irreducible states are computed by 
(1) constructing rest eigenstates of the two-body system, (2) decomposing them into irreducible 



representations under SU{2) rotations, (3) boosting the result to an arbitrary frame. The resultin 



Clebsch-Gordan coefficients in this basis are 



231] [ig] [3g] 



(Pi,A^i,P2,At2|(j,fc)p,At;^,s) 

Er/ - ^(Pl,P2) ,.-1. w, . 

''(P - Pi - P2) p N {pi,p2){l,fll,S,Hs\j,f^){jl,f^l,j2,f^2\s,f^s)X 



Yi^, {Hpi,P2))D';,^^ [RUB{p/muo), h]D';^^^ [RUB{p/muo), ^2] = (3) 
dk ^ 6{pl-pl{p,k))6{p2-p2{p,k))N{pl,p2){l,^Il,s,Hs\j,^J'){jl,fJ■'l,j2,^J'2\s,^is)x 

'J I I 



y,^,(k)D;,^^Ji?^(i?(p/mi2o),fci]I)^|^Ji?^(5(p/mi2o),A:2]. (4) 

In these expressions 

p^=p^+p^ rn?2o = V k^ = B-\plmr2^)\\{p^-p2r (5) 



ijj{px)i^{p2)K^{k) 



mi2o := 2 + m2 = 2w(A:) (7) 



and the two-body invariant mass 



is replaced by the continuous variable k := yk? . The quantum numbers / and s are kinematically 
invariant quantities that distinguish multiple copies of representations with the same mass (fc) 
and spin. For a two nucleon-system they have the same spectrum as the orbital and spin angular 
momentum operators in a partial wave representation of the nonrelativistic basis. 

An irreducible representation is also obtained by changing the order of the spin couplings in 
([3][1]) where the orbital angular momentum is first coupled to one of the spins, j2 + ^ = and then 
the result is coupled to the second spin, ji + I = j. This representation is constructed by making 
the replacements in 



'^{h,f^l,j2,f^2\s,Hs){l,fJ'l,S,fls\j,fi) 



"^{i, ^J''l, 32i ^J'2\I , iJ'iW , ^J'I , ji, \j, (8) 
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in equation or In this representation the degeneracy parameters {I, s) are replaced by (/, /). 
When we construct three-body irreducible representations by successive pairwise coupling we use 
the coupling ([3]) in the first Clebsch-Gordan coefficient and the coupling ([8]) in the second Clebsch- 
Gordan coefficient. This allows us to identify the quantum numbers of the relativistic irreducible 
basis with the quantum numbers that we used in previous nonrelativistic calculations ( 



Three-particle irreducible representations for systems of non-interacting particles can be con- 
structed by successive pairwise coupling of irreducible representations: 

(Pl, Ptl, P2, At2, P3, l^-iKJ, Q')P, A, /, J23, ^23, ^23,^23) = 

/i23 
dp23 ^ (P2,Ai2,P3,/^3|(i23,fc23)P23,/"23;^23,S23>(Pl,/^l,P23,^23|(-^,g)P,At; A,/). (9) 

M23 = -i23 

For three nucleon scattering or bound state problems it is sufficient and convenient to work in 
the three-body center of momentum frame. This simplifies the coefficients; the Wigner rotations 
in the second Clebsch-Gordan coefficient in ([9]) become the identity and the normalization factor 
N{p23,pi) 1. Both of these factors are non-trivial in the first coefficient. The form of these 



coefficients for the case that particle 1 is a spectator (for details see [20|) and using a short hand 
notation, is 

(Pl,/"'l,P2,At2)P3,At3l('^'9)P = 0,/X; A,/, j23,A:23,^23,S23> = 

5{0 - qi - q2 - q3)^r77 r 

S(.Ql -q)HH q2,q3 )-k) 

y^^(qi)y^^,(k(q2,q3)) 

Dl>^jRUB{-cii/mo23),k2{ q2,q3) [i?^(5(-qi/mo23)) ),A:3( q2,q3) )]• (10) 

where 

icii,u;iq^))=B-\P/M)p„ iG {1,2,3} (k^, ^(fci)) = ^-^-qfc/m,,)^ (11) 



^q, = M = ^^m2+q2 (12) 



i 1=1 

The important property of the states 



\ij,k)p,ti-l,s) (13) 



and 



\{J,q)P,lj.;X, I, j23 , k23 , hs , S23) (14) 
is that they transform irreducibly. The mass and spin are given by 



mi2o = 2^/m?+^ , j (15) 
M = V4m2 + 4k2 + q2 + v^m^+l? , J (16) 

respectively. 

In order to use these representations to construct dynamical representations an interaction is 
added to the two- or three-body invariant mass operator of the form 

((j, A:)p, ii---\v\---p', t,'{f, k')) = 6^^,5,,,6{p - p'){k, ■ ■ ■ \\v^ || • • , A:') (17) 

for = 2 or 

((J, q)P, fi---\V\ fi'iJ', q')) = V5jj'5(P - P')(A:, (?,••• Il^^ll • • •' , k' , q') (18) 

for = 3. Diagonalizing mi2 = mi2Q + v or M = Mq + V in the noninteracting irreducible basis 
gives simultaneous eigenstates of mi2,p, j'^ , jz for N = 2 and of M,P, J^, for = 3. In both 
the two and three-body case these eigenstates, | (ji2, Ami2)pi2, Aii2 • • • ) and \{J, Xm)P , fJ- • • • ) , ^^eie 
Xmi2 ^^'^ "^i^i eigenvalues of mi2 and M, are complete on the two and three-body Hilbert 

spaces respectively. 

The dynamical representation of the Poincare group is defined by requiring that these eigen- 
states transform like ([1]) with the mass being replaced by the mass eigenvalues Xm or Xmi2 ■ This 
representation is unitary and defines the dynamics of the system. With our choice of irreducible 
basis, {p,jcz}, the resulting irreducible representations of the Poincare group have a mass inde- 



36t | called an 



pendent representation of the three-dimensional Euclidean subgroup, which Dirac 
"instant-form dynamics" . 

For the three-nucleon case there remains the problem of how to construct realistic A^A^- 
interactions. For two-body interactions the relation 

H12 - = "2-12 = 4(k^ + m^) + AmvNN = ^mik^/m + vnn +m) (19) 

^ V ' 

h —H 

If'ni nnr 
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implies that the square of the two-body mass operator has a simple relation to the nonrelativis- 
tic rest Hamiltonian with a "realistic" NN interaction 25|] provided one identifies the spectrally 



equivalent relative momenta, k^. In the relativistic case k := B ^{p/mi2o)^{pi — P2) 

k = k(pi,p2) 

^ UJ{P2) - i^KPi) + {^{P2) - ^^{Pi)Y - V 

while in the nonrelativistic case k = S~^(p/2m)^(pi — P2)'- 

k = \ f ((Pi - i-^) - (P2 - i-m)) = kpi - P2) (21) 



2 V''^" 2m ' 2m V 2 

where i?(p/mi2o) is a rotationless Lorentz boost and Bg{Tp/2m) is the corresponding Galilean 
boost. 

The Kato-Birman invariance principle 



37 



39|] implies that the M0ller wave operators satisfy 



n±{H,Ho) := s - lim e'^*e-*^o* = s - lim e*-^(^)*e-*^(^")* = 0±(/(//), /(//q)) (22) 
where /(x) is any piecewise differentiable function of bounded variation with positive derivative 



39 1 ■ The functions f{x) = and f{x) = x^^^ satisfy the conditions of Kato-Birman theorem. 
Using equation ()22p along with the kinematic Euclidean invariance of the Hamiltonians Hy and 
Hnr gives the following relation between the two-body scattering wave operators 

0±(i/,, K,o) = n±{Hl H%) = n±{Ml M^o) = n±{Mr, Mro). (23) 

On the other hand the identification ()19p along with the reparametrization t ^ t' = Amt gives 

0±(M2,Mi) = s- lim e^^"'-^™*e-^"'-o^"** = s - lim e^^"'-*'e-'^"'-o*' = J^±„,(/i„„ /i„,o)- (24) 

Writing both wave operators as direct integrals over k^ = mhonr = {Mq — 4m^)/4 leads to the 
identifications 

n±{Hnr,Honr) = n^iKrMnr) = I ^±{\^^)d\^ = n±{Hr,Hor) = n±{Ml,M^^,) (25) 

J® 

and 

5(k2) = Ot+(k2)J^,_(k2) = J^J,,+ (k2)17„,_(k2). (26) 

The identification of the relativistic and nonrelativistic wave operators as functions of k^ ensures 
that the relativistic two-body model is fit to the same two-body 5-matrix (experimental data) as 
the nonrelativistic model provided the interactions are related by ()19p . The identification of the 
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wave operators also implies the identity of the scattering wave functions as a function of k. The 
identity of the bound state wave functions is also due to the relation p9|) . 

In our calculations we use the interaction v defined hy v := m\2 — fni20 which we construct [40I 
from the NN interaction in (jlOp by iterating 

{"1-120, v} = AmvNN - (27) 

in the irreducible plane-wave basis. Because mi2 and have the same eigenvectors, the k 
dependence of the wave functions constructed from mi2 are also identical to the corresponding 
nonrelativistic wave functions. 

The equation ()27p in the irreducible plane-wave basis has the form 



2m- 



{k,l,s\v^k',l\s') 

, I, s\v^j^^\k', I', s') ^ r „ , {k, I, s\v^\k", I", s") {k", I", s"\v^\k', I', s') 



V- f 2,, J, {k, I, s\v^\k", I", s") jk", I", s"\v^\k', I', s') 
hJ 2u.{k) + 2.{k') ^''^ 



i,{k)+u{k') ^„^„ 

The iteration converges quickly for realistic interactions [40]. While a mathematical proof 
of convergence of the iteration of (128p is lacking, the results of the iterations are easily tested 
because the resulting mi2o + v must have the same eigenfunctions as the nonrelativistic two-body 
Hamiltonian. 

We applied this approach usin g t he CD Bonn potential as the nonrelativistic interaction 



VNN{k, k'). In our previous studies 



271, 



281 ] we used the momentum transformation of ref . 24 1 



and in addition restricted to leading order terms in p/uj and v/uj expansion only 

F(k,k';,) = .(k.k') (1-^;;;^). (29) 

We checked that in most cases this simple approximation leads to practically the same results as 
the exact approach applied in the present study. 

Once the two-body mass operator is constructed, the three-body mass operator for the inter- 
acting (ij) pair is the well-defined non- linear function of the two-body mass: 



^{ij){k) = V i^ijO + Vij^ + qi + y ^2 + q|. (30) 

where Vij is embedded in the three body Hilbert space so it commutes with q^. If this is interpreted 
as the rest energy operator, the interacting pair and spectator energies are additive in the rest frame. 
This implies that the S matrix clusters properly in the rest frame - while the invariance of S in all 
frames ensures that this property extends to all inertial coordinate systems. 

11 



Pairwise interactions in the three-body system are defined by 



%)(fc) = Mf^ijm -Mo = 

{rriijo + Vij)^ + q2 - ^mf-Q + q^. (31) 

For any pair of particles these interactions commute with kinematic momentum and spin, and 
are independent on the momentum and magnetic quantum numbers. This ensures that the sum 
of the interactions has the general form p8p . 

A generalization of the method used in ()27ll28p can be used to construct y{ij){k) by iterating 



Specifically 



W^fjo + = v^ + {rriijo, v} - V(2^)(fc). (32) 



(fcr--|%)(fc)(q^)|A;',---; 



[{k, ■ ■ ■ \v'^\k', ■■■) + 2{uj{k) + uj{k')){k\v\k') 



- E /^^' • • • l%)(fc)(q')l^", • • • )k"^dk"{k'\ ■ ■ ■ |%)(fc)(q2)|A:', ■ 

111 ^11 J 



(33) 



V's" ■ 

This iteration also converges and is used to construct interactions Vij for each pair of particles. 
The three-body mass (rest energy) operator is 

M = Mo + Vi2 + V23 + V31 (34) 
where Mq is the three-body kinematic invariant mass: 



Mo = \J mj2o + q^ + Vm2 + q2 (35) 

Our relativistic Faddeev equation is based on this mass operator (|34p with two-body interactions 
constructed from the CD Bonn interaction using equations (|27p . (j28p . (|32p and (|33p . 

Finally, just like in the nonrelativistic case, there is a natural order of coupling of the irreducible 
representations for computing each pairwise interaction. The change of basis relating different 
orders of coupling is needed for the implementation of the Faddeev equation as an integral equation. 
The required basis change only changes the invariant degeneracy quantum numbers associated with 
each order of coupling. 
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((jm)P,^(a&)(c)|(/m')P',/x'((ie)(/)) = 6{P - P')<5„vA ,^^"*[(a&)(c); {de){f)] 



(36) 



The invariants R^'^[{ab){c); {de){f)] are Racah coefficients for tfie Poincare group. Tliey are con- 
structed using four Poincare Clebscli-Gordan coefficients. We compute tfiis quantity using tlie 
Balian-Brezin metfiod 



4l[ , wliere the variables associated with one order of coupling are expressed 
in terms of the variables associated with another order of coupling. The invariant coefficient, 
R^"^[{ab){c); (de)(/)], can be computed by evaluating the expression at zero momentum, averaging 
over the magnetic quantum numbers, and evaluating the resulting expression at any kinematically 



allowed set of momenta 



42 



43[. These Racah coefficients contain Wigner rotations and jacobians 



which do not appear in the nonrelativistic permutation operators. Explicit expressions are given 
in Appendix [BJ 



III. SPIN OBSERVABLES 



In a Poincare invariant quantum theory or relativistic quantum field theory the spin of a particle 
can be defined as the angular momentum that is measured in the particle's rest frame. For any 
non-zero momentum, p, different Lorentz transformations can be used to transform a particle at 
rest to frame where the particle has momentum p. Because the commutator of two rotationless 
boost generators 

[Kj,Kk] = -iejkiJj (37) 

is a rotation generator, the spin of the particle with momentum p depends on the choice of Lorentz 
transformation that transforms the particle's momentum from zero to p. In order to have an 
unambiguous definition of the spin it is necessary to choose a standard set of p dependent Lorentz 
boosts, B{p/m)^i,, that transform a particle of mass m at rest to momentum p. Then the spin of 
the particle can be unambiguously defined as the value of the spin measured in the particles rest 
frame if it is transformed to the rest frame using the standard Lorentz transformation. With this 
definition, if two spins are equal in one frame they are equal in all frames. 

The choice of standard boost is not unique because if R{p/ni) is any p dependent rotation and 

S'(p/?n) := S(p/m)i?(p/m) (38) 

then both B~^{ip/m)^y and B~^' / m)^ y both transform p to zero. 



13 



Each choice of boost leads to a different spin operator, corresponding to a different prescription 
for measuring the spin in an arbitrary frame. The rotation 

i?(p/m) = B-^{p/m)B'{p/m) (39) 



ii[22| 



that relates different boosts is called a generalized Melosh rotation 

The spin operator associated with a boost Bx{p/m) is defined as the following function of 
the Poincare generators: 

{0,3,) = B-\p/mr,W^'/m (40) 

where = ^e^^^^^paMp^ is the Pauli-Lubanski vector, Mp^ is the relativistic angular momentum 
tensor, and Bx{p/m) is the boost matrix with the parameters p/m and m replaced by the mass 
and momentum operators. 

While this quantity has the appearance of a four vector, it is not because of the operator depen- 
dence of the arguments of the boost. Instead, under Lorentz transformations the spin transforms 
like 

C/t(A,o)j,f/(A,0) = RM^,p)ix (41) 

where 

R^x{Kp) ■■= B-\Ap/m)KB^{p/m) (42) 

is the Wigner rotation associated with the boost Bx{p/m). It is a consequence of the Poincare com- 
mutation relations that the components of any of these spin observables satisfy SU{2) commutation 
relations 

[jxlijxm]— — i^lmnjxn (43) 

for any "x" . The operator is independent of the choice of boost because the generalized Melosh 
rotations leave the scalar product of two vectors unchanged. The spin defined with the textbook 
rotationless or canonical boost is called the canonical spin. 

It is natural to ask how these different types of spins are measured in the laboratory. Spins of 
isolated elementary or composite particles are measured in the laboratory through their response 
to classical electromagnetic fields. In the one-photon exchange approximation the photon cou- 
ples to matrix elements of a covariant current operator. Imposing Poincare covariance, current 
conservation, and discrete symmetries allows one to express all current matrix elements in terms 
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of an independent set of matrix elements, which have a 1-1 correspondence with invariant form 
factors. All conventional form factors can be expressed in terms of Breit frame matrix elements 
with canonical spin and a quantization axis parallel to the Breit frame momentum transfer. 

In the SL{2, C) representation canonical boosts are represented by positive Hermitian matrices. 
They have the general form 

B{p/m) = exp{a ■ p/2) (44) 

where p is the rapidity of the Lorentz transformation and cr are the Pauli matrices. In this paper all 
of our spins are canonical spins. The 50(1, 3) representation of canonical boosts are the standard 
rotationless boosts. 

Given a definition of the form factors in terms of independent current matrix elements in a given 
basis, it is also possible to express them in terms of current matrix elements in any other standard 
frame using any other basis [45]. For example the expression in terms of Breit frame canonical 
spin matrix elements, can be replaced by a different independent set of laboratory frame helicity 
spin matrix elements. In quantum field theory the choice of boost is built into conventions used 
to define the Dirac spinors. The relation of the invariant form factors to current matrix elements 
with different choices of spin determines the relationship between different spin observables and 
experiment. 

The spin degrees of freedom of the asymptotic incoming or outgoing particles are most con- 
veniently expressed in terms of traces of density matrices, which is a reflection of the fact that 
realistic initial and/or final states are generally not pure states. Scattering spin observables in 
cross sections are formally defined by [3J] 

^ Tr(5,r5.Tt) 
^ ' Tt{TT^) ^ ^ 

where T is the invariant scattering amplitude for the reaction under consideration. The connection 

between the invariant scattering amplitude defined in the particle data book 46|] and the transition 



amplitudes constructed by solving our formulation of the relativistic Faddeev equations is given in 
The quantities S have the form 

S = ^SaSa (46) 

where the index runs over all Ni or Nf initial or final sets of magnetic quantum numbers, Sa are 
a basis for Ni x Ni or Nf x Nj matrices that are orthonormal with respect to the trace norm and 



Sn are constant coefficients 



3i . 
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If the initial and final asymptotic states are represented in a canonical spin basis, then the 
magnetic quantum numbers that appear in the invariant amplitudes T transform with Wigner 
rotations under Lorentz transformations. The result is that the spin observable (O) will not be 
invariant unless the matrices Sa or coefficients Sa are defined to transform in a manner that leaves 
the observable invariant. 

Any spin observable can be made Lorentz invariant, by defining the invariant observable as its 
value in a given frame if it is transformed to the frame with a specific Lorentz boost. This can be 
used to get an invariant definition of the vector or tensor polarizations. 

In this paper invariant spin observables are defined to be the values of the observable in the 
laboratory frame (rest frame of the target). To evaluate the corresponding spin observable it is 
only necessary to evaluate the expression 

^ Tr(g,rg.Tt) 

^ ' Tr(rrt) ^^'^ 

for the values of the invariant amplitudes with laboratory kinematics. 

This observable is equal its value evaluated in other frames using the formula 

Tr{SfD^MD'S,D'^M^D) 
Tr(MMt) 

where the invariant amplitudes are evaluated in the other frame and D and D' are products of 
Wigner D-functions of the Wigner rotations associated with the boost from the lab frame to the 
other frame. 

Our specific interest in this paper is the observable Ay with polarized incoming nucleon. The 
convention used to define Ay is the Madison convention, where the laboratory frame scattering 
plane is in the xz plane. The observable Ay, is defined as 

_ TV(TK X /.)rt) 
{^y) - Tr(rrt) ^ ^ 

Because the Wigner rotation for canonical boost along the direction of a particle's momentum is 
the identity, ay is unchanged so Ay can also be evaluated in the cm. frame without making any 
compensating Wigner rotations. 

IV. FADDEEV EQUATION 

The nucleon-deuteron scattering with neutron and protons interacting through a NN interaction 
vnn alone is described in terms of a breakup operator T satisfying the Faddeev-type integral 
equation 



341, 



48] 
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T\<P) = tP\cl))+tPGoT\cl)). (50) 

The two-nucleon (2N) t-matrix t results from solving the Lippmann-Schwinger equation with 
the interaction vnn- The permutation operator P = P12P23 + ^i3-F23 is given in terms of the 
transposition Pij which interchanges nucleons i and j. The incoming state \(j)) = \(io)\(f)d) describes 
the free nucleon-deuteron motion with relative momentum qo and the deuteron state vector \(j)d)- 
Finally Go is resolvent of the three-body center of mass kinetic energy. 

The elastic nd scattering transition operator U is given in terms of T by {34, 



48| 



U = PGq^+PT. (51) 

This is our standard nonrelativistic formulation, which is equivalent to the nonrelativistic 3N 
Schrodinger equation plus boundary conditions. The formal structure of these equations in the 
relativistic case remains the same but the ingredients change. As explained in 2^ the relativistic 



3N rest Hamiltonian (mass operator) has the same form as the nonrelativistic one, only the mo- 
mentum dependence of the kinetic energy changes and the relation of the pair interactions in the 
three-body problem to the pair interactions in the two-body problem changes. Consequently all 
the formal steps leading to Eqs. (|50p and ()5ip remain the same. 

The free relativistic invariant mass of three identical nucleons in their cm. system has the form 
(j35p while the free two-body mass operator has the form ([7]). 

As introduced in [l9(] the pair forces in the relativistic 3N 2-1-1 mass operator are given by 
()3ip where V = V{q^) reduces to the interaction v for q = 0. 

The transition matrix that appears in the kernel of the Faddeev equation (jSOp is obtained by 
solving the "Lippmann-Schwinger equation" , which must be solved as a function of q^ 



^(2u;(k ') 2 + q 2 - V(2t^(k ") 2 + q 2 + 
The input two-body interactions are computed by solving equation ([32]) and ([33]) . 

The new relativistic ingredients in Eqs. (|50p and (15ip will therefore be the i-operator (I52[) (ex- 
pressed in partial waves) and the resolvent of the 3N invariant mass 

where Mq is given by Eq. (I35p and E is the total 3N cm. energy expressed in terms of the initial 
neutron momentum qo relative to the deuteron 



E = ^ (Mrf)2 + qg2 + ^^2 + q^2^ (54) 
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and Md is the deuteron rest mass. 

Currently Eq. (|5U|) in its nonrelativistic form is numerically solved for any NN interaction using a 



This turns Eq. (j50j) 
20( 1 ■ in the relativistic 



momentum space partial-wave decomposition. Details are presented in ref. [45 
into a coupled set of two-dimensional integral equations. As shown in 
case we can keep the same formal structure, though the permutation operators are replaced by 
the corresponding Racah coefficients ()36p for the Poincare group. These coefficients include both 
jacobians and Wigner rotations that do not appear in the nonrelativistic permutation operators 



35|. These coefficients are computed in the Appendix [Bl using methods that we have applied to 
compute the nonrelativistic permutation operators. 

In the nonrelativistic case the partial wave projected momentum space basis is 

\pq{ls)3[x\)IJ{t\)T), (55) 



where p and q are the magnitudes of standard Jacobi momenta (see 3J, l35|), obtained by trans- 
forming single particle momenta to the rest frame of a two or three-body system using Galilean 
boosts, and {ls)j two-body quantum numbers with obvious meaning, (A1/2)I refer to the third 
nucleon (described by the momentum q), J is the total 3N angular momentum and the rest are 
isospin quantum numbers. In the relativistic case this basis is replaced by the irreducible plane 
wave states defined in ([10]). 

The basis states (jlOp are used for the evaluation of the partial wave representation of the 
permutation operator P with Wigner rotations of spin states for nucleons 2 and 3 included. In the 
relativistic case we adopt the following short-hand notation for the irreducible three-body states, 
which also includes isospin quantum numbers coupled in the same order: 

\k,q,a) := \kq{ls)j{\]^)IJ{t]^)T) = |( J, g)P, /i; A, /, J23, A:23, ^23, S23)|(t^)T) (56) 

Equipped with that, projecting Eq. (|50|) onto the basis states \k,q,a) one encounters like in the 



nonrelativistic notation 



i{kqa\P\k' q a')i = i{kqa\k'q'a')2 +1 {kqa\k'q'a')3 = 2 i{kqa\k'q'a)2- (57) 

This is evaluated by inserting the complete basis of states |pi, ^ui, p2, /U2, P3, /^s > and using . It 
can be expressed in a form which resembles closely the one appearing in the nonrelativistic regime 



i{k q a\ P \k' q' a)i = J dx 



1 ^ 6{k-Tri) 5{k' -^2) 



A:2 A:'2 
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where all ingredients are given in the Appendix [Bj 

Due to short-range nature of the NN interaction it can be considered negligible beyond a certain 
value jraax of the total angular momentum in the two-nucleon subsystem. Generally with increasing 
energy jmax will also increase. For j > jmax we put the t-matrix to be zero, which yields a finite 
number of coupled channels for each total angular momentum J and total parity vr = (— )'^'^ of 
the 3N system. To achieve converged results at our energies we used all partial wave states with 
total angular momenta of the 2N subsystem up to jmax = 5 and took into account all total angular 
momenta of the 3N system up to J = 25/2. This leads to a system of up to 143 coupled integral 
equations in two continuous variables for a given J and parity. 



V. RESULTS AND DISCUSSION 



The subject of the present study is to investigate the influence of relativity on the nd elastic 
scattering nucleon analyzing power Ay at low energies. We define the invariant observable Ay to 
be the value measured in the laboratory frame (target at rest). 

To this aim we solved Faddeev equations at a number of the incoming neutron lab energies £"„ = 
5 MeV, 8.5 MeV, and 13 MeV. To check the energy dependence of the effect we added two additional 
energies En = 35 MeV and 65 MeV. In order to see the importance of specific relativistic features 
we solved the equation in the relativistic case with and without Wigner rotations. This allowed us 
to see which effects, dynamical corrections, which are induced by the momentum dependence of the 
two-body force together with kinematical relativistic effects coming from the use of the Poincare 
Jacobi variables, or Wigner rotations, dominate for Ay. 

Figs H] and [2] illustrate the results. When only dynamical effects are taken into account (Wigner 
rotations neglected) then at low energies of the incoming neutron the relativistic and nonrelativistic 
predictions are practically the same with the exception of the angular region close to the maximum 
of Ay where the relativistic prediction is ~ 2% below the nonrelativistic one (see Fig(T]). This small 
effect disappears at higher energies (see Figl2]). Including Wigner rotations lowers significantly 
the values of Ay in a large region of angles around the maximum. The changes in the maximum 
are up to ~ 10%. Again, when the energy of the incoming neutron increases nonrelativistic and 
relativistic predictions are practically identical. 



The large changes of Ay occur in a region of energies where t 



lis observable is extremely sensitive 



to changes in ^Pq, ^Pi, and F2 NN force components [49|. At energies where this sensitivity 
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dies out also the relativistic effects for Ay become negligible. This allows to conclude that the large 
effects seen for Ay at low energies are due to amplification of changes of the ^Pj contributions due 
to relativity by a large sensitivity of Ay to P- waves. 

In a recent study [29] the changes of Ay due to relativity by 10% at E^^^ = 3 MeV have been 
reported. Very probably the opposite sign of the effect found in that study can be attributed to 
the impulse approximation used when calculating Ay. 



VI. SUMMARY AND OUTLOOK 



We numerically solved the 3N Faddeev equation for nd scattering including relativistic kine- 
matics, dynamical relativistic effects and Wigner rotations at the neutron lab energies i?^"* = 5, 
8.5, 13, 35, and 65 MeV. As dynamical input we took the nonrelativistic NN potential CD Bonn 
and generated in the 2N cm. system an exactly on-shell equivalent relativistic interaction v, by 
solving numerically nonlinear quadratic equation relating matrix elements of the nonrelativistic 
and relativistic potentials. We checked that the approximate procedure using an analytical scale 
transformation of momenta applied in our previous studies provides practically the same results 
as the present exact approach. In addition a similar nonlinear equation (Eq. (j32p ) was used to 
generate the momentum dependent two-body interaction embedded in the 3-particle Hilbert space. 

We found that at low energies the effects of Wigner rotations are most important for the 
analyzing power. They lower the maximum of Ay by up to ~ 10%. The dynamical relativistic 
effects are of minor importance for Ay and provide small changes of Ay in a region close to its 
maximum. They lower Ay by only ~ 2%. The relativistic effects disappear at higher energies. 

Wigner rotations are negligible for the cross section and all other spin observables in elastic 
Nd scattering with exception of four low energy spin correlations Cz,yz, Cx,xy, Cy^xx-yy, Cy^yy and 
four low energy spin transfer coefficients: from deuteron to deuteron Kz ^ and Kx ^ , and from 
deuteron to neutron Kyz and Kxy Very probably this can be traced back as in the case of Ay to 
their sensitivity to '^P-waves. 

These results shed new light on the low-energy analyzing power puzzle. It is known that the 
existing discrepancies between Ay data and theoretical predictions based on NN potentials only 
cannot be removed when current three-nucleon force, mostly of 27r-exchange character [l^, |ll| are 
included in the nuclear Hamiltonian. This indicated that additional 3N forces should be added 
to the 27r-exchange type forces. Such forces provided by yPT in NNLO and NNNLO orders are 
expected to provide the solution for the Ay puzzle [sj] • It seems that in view of the present 
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FIG. 1: (color online) The nucleon analyzing power Ay for nd elastic scattering at various lab energies E^^^ 
of the incoming neutron. The dotted line is the result of the nonrelativistic Faddeev calculation with the 
CD Bonn potential. The relativistic predictions without and with Wigner spin rotations are shown by the 
dashed and solid lines, respectively. The nd data at 5 MeV and 8.5 MeV are from ref. [s^ and at 13 MeV 
are from ref. 15 il l. 
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FIG. 2: (color online) The nucleon analyzing power Ay for nd elastic scattering at £"„ = 35 MeV and 
65 MeV. For description of lines see Fig[T] All theoretical predictions are practically overlapping. The pd 



data at 35 MeV are from ref. 



53 and at 65 MeV from [53 1 



result they must increase the maximum of Ay stronger than expected up to now. However, Ay is 
a very sensitive observable and our approach, using the irreducible {p,jcz} basis, is only one of 
many possible basis choices that lead to Poincare invariant dynamical theories that are two-body 
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scattering equivalent, which may give different three-body predictions. 
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APPENDIX A: SPINORS 



In our calculations SL(2, C) matrices are used to represent Wigner rotations and Lorentz trans- 
formations. In this way spin algebra is reduced to working with 2x2 complex matrices, and there 
is no need for representations in terms of Euler angles. The relevant relations are give below. 

The coordinates of a four-vector can be labeled by the 2x2 Hermitian matrix P: 

(„0 _|_ r)3 pi _ \ 
1 2 3 ' 
+ ip^ p" — p"^ I 

where are the identity and the three Pauli matrices. The components of p^ can be extracted 
from the matrix P using: 

p'^ = ^Tvia.P) = ^Tr{Pa,). (A2) 

Because 

det(P) = - {p f = -ri^^p'^p'' = (A3) 

and 

P = (A4) 

for realp'^, it follows that any linear transformation that preserves the Hermiticity and determinant 
of P is a real Lorentz transformation. It is easy to show that if A is a complex 2x2 matrix with 
det(A) = 1 then the transformation 

P^P' = APA^ (A5) 
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has both of these properties: 

det(P') = det(P) and P = ^ P' = P'^ . (A6) 

The most general 2x2 matrix with determinant 1 can be written as 

A = ±e^^-^ (A7) 

where z = 6 — ip. If p = then A = U{6) is unitary and corresponds to an SU{2) rotation through 
an angle \6\ about the axis. If = then A is a positive Hermitian matrix that corresponds to 
a rotationless (canonical) Lorentz boost with rapidity \p\ in the direction p. 
The rotation U{9) is given by 

6 9 

A^U{e) = aoCOs{-)+icr-esm{-). (A8) 
The axis of rotation can be extracted using 

= -i ^("^W) (A9) 

and the angle of rotation can be extracted from 

The rotationless boost that transforms a particle of mass m at rest to total momentum p can 
be labeled by the final four velocity q := p/m: 

A ^ Bid) = aocosh(^) + a- ■ psinh(|) (All) 

where p = pp is the rapidity of the Lorentz boost which is related to q by 



and 



B{q) satisfies 



-^^(f) = VV = V^ (A12) 



sinh(J) = , '"^^ = , (A13) 



5(q)| U^(q)= \ \\ \\=P (A14) 

0ml \ P + ip P — P 
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where 

p° = uj{p) = A/m2 + p • p. (A15) 

The rotationless boost above is called the canonical boost. The inverse transformation is obtained 
by reversing the sign of q or p. 

The Wigner D functions are homogeneous polynomials of degree 2j in the coefficients of the 
5L(2,C) matrices Aij-. 

D'^iA, := f: l(i + - - ^j^>-M;,A;r'-"-4r- (A16) 

(j + /X — a)!a!(a — /i + — 1/ — aj! 

These are representations of both SL{2,C) and SU{2). The j = 1/2 representation is just the 
matrix A. 

APPENDIX B: PERMUTATION OPERATOR 

Using Eq. (fTOj) twice for the bra state a| and the ket state \k\q',a')2 one gets for the 

matrix element of the permutation operator in our partial wave basis: 

i{k ,q ,a\ P \k' ,q,' a')i = 2 i{k ,q ,a \ k' , q' , a')2 = 

2 E E E E E 

mim^m^ ^42/^3/^8 /i;MA/i//i /ij/^sMs' Mj'A'a'Mi'M' 

(Ama^,"ii, {j,fi,i,fii\j,M) {^,^J.2,^,^i3\s,^^s) {i,ni,s,^is\j,fi) 

{>^',^J'\',^,'m2,\I',^J'I') {j',n',l',m'\J,M) {^,fJ.2,^,l^3\s',fls') {l',IJ'l',s',IIs'\j',n') 

dq fiq' 



iV(q',-q-q') Ar(-q-q',q) 
<5(fc-|k(q',-q-q')| ) 5( A:'-|k( -q-q',q 



with 



k2 k'2 

y^/(q) y/;(k( q ',-q - q ') ) Y^^iq ') y^;,(k( - q - q ',q ) ) 
I?i%[i?«,(5(-q/2cum(k)),(k( q',-q-q') ,c^„(k))] 
£'i3M3[^«.(^(-q/2c^m(k)),(-k( q',-q-q' ) ,a;„(k))] 
D|^^,[ii^(i3(-q' /2c^„(k)),(k( -q-q',q) ,a;„(k))] 

Dl^^,[R^{B{-q' /2coM),{-H - q - q q ) , Wm(k))] 

1 < ( (4)^ I >2, (Bl) 



k(q',-q-q') = <i ' + + Viiq, q' , x)) 
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k( - q - q ',q 



q- 'i'i- + y2iq,q',x)). 



(B2) 



In Eq.(lB2]) occur 



Eq' — E, 



q+q' 



^q' + ^q+q' + ^(^q' + ^q+q')^ " q ^ 

with x = q-q', y2{q, q', x) =yiW, q, x) and E^ = ujm{<l )• 



(B3) 



Proceeding as in 



41 



42 



43] one gets the following expression for the matrix element of the 



permutation operator P: 

i{k ,q a\ P \k' , q' a')i 

with 



\ 5{k- TTi) 6{k' - TT2) 

dx- 



-1 



Ni{q,q>,x) N2iq,q',x) 



Gaa'iq^ q'^ x) 



(B4) 



TTi = \lq'^ + -q^{l + yi)^ + qq'x{l+yi) 



= yg^ + + 2/2)^ + g'g''a;(l + ^2) 

Ni{q,q',x) = iV(q',-q-q') 
N2{q,q\x) = iV(-q-q',q), 



(B5) 



and 



Gaa'iqq'x) = ^ ^ {-If Stt'SmtMj,, V it 



4vr3/2 



-1/2 1/2.^ 



1/2 T t' 



^2" 2' 



X](2''^2' 2'^3N''/^2 + ^3) 



^(AO^, mi, |/, mi)Z);/f^, [i2^(S(-q72a;^(k')), (-k', w„,(k'))] 

mi 

^{l, /i - //2 - At3, ^2 + n){-)''~^^~^^Yi_(^f,_^^_^.^{p){j, /i, I, mi I J, /i + mi) 
^(A',;U + mi - /i' - m2, -,m2|l',/i + mi - /x') 

+mi— /i'— m2Vq J- 



(B6) 
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We use standard notation I = 21 + 1. It is assumed that z-axis is along q and the momentum q ' 
Ues in the x-z plane what leads to the following components of the q, q k, and k ' vectors 



Though t 
used like in 



q = [0,0, g] 

q' = [q'\/ 1 — x^, 0, q'x] 
k = [gVl - x^, 0, q'x + ^q{l + yi{q, q' , x))] 

k' = [-^q'il + y2iq, q', x))Vl - x\ 0, -q - ]^q'{l + y2{q, q' , x))x] (B7) 

re direct evaluation of Euler angles as arguments of the Wigner D-functions could be 
201 1 ■ here we use the SL{2,C) representations of Lorentz transfomations discussed in 



Appendix O This leads to 

D^/^[R.UB{-ci/Mo),{k,m))] = S(p2/m)S(-q/Mo)S(k/ 



where 



Eo + Mq / u}{k) + m 
2Mo \jio{p2)+m 
k ■ q 



^2Mo{Eq + Mo)(tj(fc) + m){Lo{p2) + m) 



+ ik X q • cr 



^2Mq{Eq + MQ){uj{k) + m){uj{p2) + m) 



(B8) 



Eq 

Mo 

P2(k, -q) 



M2 + g2 
2Lo{k) 



kq 



Mo(^o + Mo)' 



(B9) 
(BIO) 
(Bll) 
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